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ENUMERATION OF A DUAL SET OF STIRLING PERMUTATIONS BY 

THEIR ALTERNATING RUNS 

SHI-MEI MA AND HAI-NA WANG 


Abstract. In this paper, we count a dual set of Stirling permutations by the number of alter¬ 
nating runs. Properties of the generating functions, including recurrence relations, grammatical 
interpretations and convolution formulas are studied. 
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1. Introduction 


Let [n] = {1, 2,... , n}. The Stirling number of the second kind is the number of ways to 

partition [n] into k blocks. Denote by D the differential operator and let id = xD. It is well 

known that 


k =1 J 


z k D k . 


Let 


r(x) = 


VTT 


X 


-y/l — x 

By induction, one can easily verify that there are positive integers T(n,k), k € [2n — 1], such 


that 


Y?kJi T (n, k)x k 


$ n (r(x)) = 


for n > 1. 


(1 — x) n (l+x) n 1 v / l — x 2 
It is clear that the numbers T(n,k ) satisfy the initial conditions T(l,l) = 1 and T(l,k) = 0 
for k / 1. Let T n (x ) = f 1 T{n,k)x k . Using $ n+1 (r(x)) = i9('i? n (r(x))), we get that the 
polynomials T n (x ) satisfy the recurrence relation 

T n+ i(x) = (2 nx + 1 )xT n (x) + x(l - x 2 )T^(x) (1) 

for n > 0, with the initial values Tq(x) = 1 and T\(x) = x. In particular, 

T n { 1) = -T n+1 (-1) = (2 n - 1)!! for n > 1. 

Equating the coefficients of x k on both sides of ([Tj) , we get that the numbers T(n, k ) satisfy the 
recurrence relation 


T(n + 1, k) = kT(n, k ) + T(n, k — 1) + (2 n — k + 2 )T(n, k — 2). (2) 

The motivating goal of this paper is to find a combinatorial interpretation of the numbers T(n, k). 

Let & n denote the symmetric group of all permutations of [n], where [n] = {1,2,... ,n}. Let 
7 r = 7 t(1)7t( 2) • • • 7r(n) € & n . An interior peak in 7r is an index i € {2, 3, ..., n — 1} such that 
7r(i—1) <7r(z) > 7r(z+l). A left peak in ir is an index i € [n—1] such that 7r(i—1) < 7r(i) > 7r(i+l), 
where we take 7r(0) = 0. Let ipk(7r) (resp. lpk(7r)) be the number of interior peaks (resp. left 
peaks) in it. We say that 7r changes direction at position i if either 7r(i — 1) < tt(i) > 7r (i + 1), 
or 7r(i — 1) > ir(i) < 7r(i + 1), where i € {2,3,..., n — 1}. We say that 7r has k alternating runs if 
there are k — 1 indices i such that 7r changes direction at these positions. Denote by altrun (7r) 
the number of alternating runs in 7r. 


l 
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Define 

W n {x) = Y z ipkW > W n (x) = Y z lpk(7r) , Rn{x) = Y x altrun(7r) . 

irG&n 7rG6 n 7re6 n 

From m Corollary 2, Theorem 3], we get 

(1 2aT )2 ^ n(x) = + ^n(s 2 )- 

Let Rn(x) = Ylk=\ R( n i k)x k . The study of alternating runs of permutations was initiated by 
Andre [2], and he proved that the numbers R(n , k) satisfy the recurrence relation 

R(n, k ) = kR(n — 1, k) + 2R(n — l,k — 1) + (n — k)R(n — 1, k — 2) (3) 

for n,k > 1, where i?(l, 0) = 1 and 12(1, A;) = 0 for k > 1. It follows from Q that the polynomials 
R n {x) satisfy the recurrence relation 

R n + 2 (x) = x(nx + 2)R n+1 (x) + x (l - x 2 ) R' n+1 (x), (4) 

with the initial value R\(x) = 1. Recall that a descent of a permutation it € & n is a position i 
such that 7r(i) > ir(i + 1). Denote by des (7r) the number of descents of it. Then the equations 

A n (x)= Y x des ^ +1 = Y(l)x\ 

it e&n k= 1 ' ’ 


define the Eulerian polynomial A n {x) and the Eulerian number The polynomial R n (x) is 
closely related to A n (x): 


Rn{x) = 


1 + x 


71—1 


(1 + w) n+l A 


1 — w 
1 + w 


w = 


1 — X 

1 + x' 


( 5 ) 


which was first established by David and Barton [71 157-162] and then stated more concisely 
by Knuth El p. 605]. There is a large literature devoted to the polynomials R n {x) (see [221 
A059427]). The reader is referred to [Hill] for recent results on this subject. 

In [5], Carlitz introduced C n (x) defined by 



C n (x) 

(1 — x) 2k+1 ’ 


and asked for a combinatorial interpretation of C n {x). Riordan m noted that C n (x) is the enu¬ 
merator of trapezoidal words with n elements by number of distinct elements, where trapezoidal 
words are such that the i-th element takes the values 1,2,..., 2* — 1. Gessel and Stanley [8j gave 
another combinatorial interpretation of C n (x) in terms of descents of Stirling permutations. A 
Stirling permutation of order n is a permutation a = a(l)cr(2) • • • a(2n — l)cr(2n) of the multiset 
{1,1, 2, 2,... ,n, n} such that for each i, 1 < i < n, all entries between the two occurrencies of 
i are larger than i. Denote by Q n the set of Stirling permutation of order n. For a € Q n , let 
ct(0) = a(2n + 1) = 0, and let 


des (a) = #{* | a(i) > cr(i + 1)}, 
asc (a) = #{* | a(i — 1) < 
plat (a) = #{* | a(i) = a(i + 1)} 


denote the number of descents, ascents and plateaux of a, respectively. Gessel and Stanley [8] 
proved that 

C n (x) = Y xAeSa - 

O-eQn 
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Bona [31 Theorem 1] introduced the plateau statistic on Q ni and proved that descents, ascents 
and plateaux are equidistributed over Q n . The reader is referred to mmm for recent progress 
on the study of statistics on Stirling permutations. 

In the next section, we show that T n (x) is the enumerator of a dual set of Stirling permutations 
of order n by number of alternating runs. 

2. Combinatorial interpretation of T(n,k) 

Let a = cr(l)(j(2) • • • er(2n) € Q n ■ Let <f> be the injection which maps each first occurrence of 
entry j in a to 2 j and the second j to 2 j — 1, where j € [n]. For example, <£(221331) = 432651. 
The dual set <f>(Q n ) of Q n is dehned by 

${Qn) = {7T | a £ Qn, <£(<t) = 7T } . 

Clearly, <f>(Q n ) is a subset of © 2 n- For n € <£(<2 n ), the entry 2 j is to the left of 2j — 1, and all 
entries in ir between 2 j and 2 j — 1 are larger than 2j, where 1 < j < n. Let ab be an ascent in 
a, so a < b. Using <f>, we see that ab maps into (2a — 1)(26 — 1), (2a — 1)(26), (2a)(26 — 1) or 
(2a)(2b), and vice versa. Note that asc (<r) = asc (<£(<r)) = asc(7r). Therefore, we have 

c n ( x ) = * asc(7r) - 

7T6'I’(Qn) 

It should be noted that ir € <£(Q ?l ) always ends with a descending run. We now present the 
following result. 

Theorem 1. We have 


T(n , k) = #{ 7 r € $(Q n ) | altrun (t) = k}. 

Proof. There are three ways in which a permutation ir € <£(Q n + 1 ) with altrun (-k) = k can be 
obtained from a permutation a € 4 > (Q n ) by inserting the pair (2n+2)(2n+l) into consecutive 
positions. 

(a) If altrun (a) = k, then we can insert the pair (2n + 2) (2 n + 1) right before the beginning 
of each descending run, and right after the end of each ascending run. This accounts for 
kT(n,k) possibilities. 

(b) If altrun (a) = k — 1, then we distinguish two cases: when a starts in an ascending run, 
we insert the pair (2n + 2)(2n + 1) to the front of a; when a starts in an descending run, 
we insert the pair (2 n + 2)(2n + 1) right after the first entry of a. This gives T(n, k — 1) 
possibilities. 

(c) If altrun (a) = k — 2, then we can insert the pair (2 n + 2)(2n + 1) into the remaining 
(2n + l) — (k— 2) — 1 = 2n — k+2 positions. This gives (2n — k+2)T(n,k—2) possibilities. 

Therefore, the numbers T(n, k) satisfy the recurrence relation Q, and this completes the proof. 

□ 


Define 

M n (x) = ® ipkW : N n (x) = a;lpkW - 

7re4>(Sn) 7rg$(Q„) 

It follows from m Theorem 4] that M n (x) = x n N n (b). Moreover, from [16, Theorem 5], we 
have 

(1 + x)T n (x) = xM n (x 2 ) + N n (x 2 ). 

We now recall some properties of N n (x). Let N n (x) = N(n, k)x k . Apart from counting 

permutations in the set <£(Q n ) with k left peaks, the number N(n,k) also has the following 
combinatorial interpretations: 
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(mi) Let e = (ei, e^-, ■ ■ ■, e n ) € Z n , and let I n k = {e € Z n |0 < e* < (i — l)k}, which known as 
the set of n-dimensional k-inversion sequences (see (20]). The number of ascents of e is 
defined by 


asc (e) 


# 


/i:l<i<n—1 


Cj &i +1 1 

(i — 1 )k + 1 ik + 1 J 


Savage and Viswanathan ISP discovered that N(n, k) = #{e € / nj 2 : asc (e) = n — A:}. 
(m 2 ) We say that an index i G [2n—1] is an ascent plateau of n € Q™ if ir(i— 1) < tt( i) =7r(z+l). 
The number ZV(n, A:) counts Stirling permutations in Q n with k ascent plateaux (see [161 
Theorem 3]). 

(m 3 ) The number N(n,k ) counts perfect matching on [2n] with the restriction that only k 
matching pairs with odd minimal elements (see EE]). 

The polynomials ZV n (x) satisfy the recurrence relation 


N n+ i(x ) = (2n + 1 )xN n (x) + 2x(l — x)N' n (x ) 


with initial value ZVo(x) = 1. The first few of N n (x ) are 

ZVi(x) = x, N 2 (x) = 2x + x 2 , N$(x) = 4x + 10 x 2 + x 3 , N±(x) = 8 x + 60x 2 + 36x 3 + x 4 . 

The exponential generating function for N n (x ) is given as follows (see m Section 5]): 

N (x, Z ) = ^2 N n(x)^ = ^Y~ x ~ 2 ^x- x y ( 6 ) 


A polynomial /(x) = Y^k=Q a k xk is symmetric if = a n _fc for all 0 < k < n, while it is 
unimodal if there exists an index 0 < m < n, such that 


0*0 ^ ^ ^ ®m-1 A dm ®m+l ^ ^ ®n- 

Theorem 2. T/ie polynomial T n (x) is symmetric and unimodal. 

Proof. It is immediate from ([2]) that T n {x) is a symmetric polynomial. We show the unimodality 
by induction on n. Note that Ti(x) = x,T 2 (x) = x + x 2 + x 3 and T 3 (x) = x + 3x 2 + 7x 3 + 3x 4 + x 5 
are all unimodal. Thus it suffices to consider the case n > 3. Assume that T n (x) is symmetric 
and unimodal. For 1 < A; < n + 1, it follows from (] 2 ]) that 

T(n + 1, A:) — T(n + 1, k — 1) = (A; — l)(T(n, k ) — T(n, k — 1)) + ( T{n , k — 1) — T{n, k — 2)) 

+ (2n — k + 2 )(T(n, k — 2) — T(n, k — 3)) + (T(ra, A:) — T(n , k — 3)) 
> 0 , 

where the inequalities are follow from the induction hypothesis. This completes the proof. □ 


In the next section, we present a grammatical interpretation of T n (x). 


3. Grammatical interpretations 

The grammatical method was introduced by Chen [ 6 ] in the study of exponential structures in 
combinatorics. For an alphabet A, let Q[[A]] be the rational commutative ring of formal power 
series in monomials formed from letters in A. A context-free grammar over A is a function 
G : A —» Q[[A]] that replace a letter in A by a formal function over A. The formal derivative 
ZD is a linear operator defined with respect to a context-free grammar G. More precisely, the 
derivative D = Dq. Q[[A]] —>• Q[[A]] is defined as follows: for x € A, we have D{x) = G(x)\ 
for a monomial u in Q[[A]], D{u ) is defined so that ZD is a derivation, and for a general element 
q € Q[[A]], D(q) is defined by linearity. 
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The hyperoctahedral group B n is the group of signed permutations of the set ±[n] such that 
7r(—z) = — 7r(z) for all i, where ±[n] = {±1, ±2,..., ±n}. For each n G B n , we define 

desA(vr) := #{i € {1,2,... ,n - l}|7r(i) > tt( i + 1)}, 
des b(7t) := #{i G {0,1,2 ,n- l}|vr(i) > it (i + 1)}, 

where 7r(0) = 0. Following [1], the flag descent number of ir is defined by 

f 2des a(vt) + 1, if 7r(l) < 0; 
fdes (7TJ := < 

l2desA(7r), otherwise. 

Let 


n 


B n {x) 

— ^.des B (tt) 

= ^2B(n,k)x k , 


ir£B„ 

/c=0 



2 n 

S n (x) 

— x fdes ( 7r ) = 

= 22 Sin, k)x k ~ 1 


7T £B n 

k= 1 


The polynomial B n (x ) is called an Eulerian polynomial of type B, while B(n,k ) is called an 
Eulerian number of type B (see [221 A060187]). It follows from [I] Theorem 4.3] that the 
numbers S(n,k ) satisfy the recurrence relation 

5(n, k) = kS(n — 1, k) + 5(n — 1, k — 1) + ( 2 n — k + 1 )S(n — 1, k — 2) 

for n, k > 1, where 5(1,1) = 5(1,2) = 1 and 5(1, k) = 0 for k > 3. The polynomials S n (x) is 
closely related to the Eulerian polynomial A n (x ): 

5 n (ic) = -(1 + x) n A n (x) for n > 1, 
x 

which was established by Adin, Brenti and Roichman [11. It should be noted that S n (x ) and 
A n (x) are both symmetric. 

Consider the context-free grammar 

A = {x, y, z }, G = {x p(x, y, z),y q{x, y, z),z -A r(x, y, z)}, 

where p(x,y, z),q(x,y, z) and r(x,y,z ) are polynomials in x,y and z. The diamond product of 
z with the grammar G is defined by 

Go 2 = {ig p(x, y, z)z, y ->• q(x, y, z)z, z ->• r(x, y, z)z}. 

We now recall two results on context-free grammars. 

Proposition 3 (p31 Theorem 6 ]). If 

G = {x ->• xy,y -A- yz,z y 2 }, (7) 


then 

n 

D n (x 2 ) = x 2 ^ R{n + 1, k)y k z n ~ k . 
k =0 

Setting x = z = l, we have D n (x 2 )\ x=z= \ = R n+ \{y). 

Proposition 4 m Theorem 10]). Consider the context-free grammar 

G' = {x xyz , y -A- yz 2 , 2; -A y 2 z}, 


( 8 ) 
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which is the diamond product of z with the grammar G defined by 0 - For n > 1, we have 


2 n 


D n (xy)=xY,S(n,k)y 2n - k+1 z\ 
k =1 
n 

D n (yz ) = Y J B{n 1 k)y 2n ~ 2k+1 z 2k+ \ 

k =0 
n 

D n {y) = ^2N(n,k)y 2n ~ 2k+1 z 2k , 
k =1 
n 

D n (z) = J2 N{n, n-k + 1 ) y ^- 2 k +2 z 2 k-i^ 


k =1 


D n (y 2 ) = 2 n Y J 


k =1 


2n—2/c+2^,2/c 


We can now conclude the following result. 

Theorem 5. Let G' be the context-free grammar given by (| 8 j). Then for n > 1, we have 


2n—l 


Jln—k 


D n {x ) =x^2 T(n,k)y k z 

k =1 

D n (x 2 ) = 2x 2 (y + zr- 1 f2(l)y‘ 

k =1 ' ' 


k ^n—k+1 


Setting x = z = 1, we /lave L> n (x)| a;= ^ = i = T„(y) and D n (x 2 )\ x=z =i = 2(1 + y) n 1 A n (y). 

Proof. Note that D(x) = xyz and D 2 (x) = xyz 3 + xy 2 z 2 + xy 3 z. For n > 1, we dehne t(n, k ) by 

£> n (s) =x^t(n,/t)y fc 2 2n - fc . 


fc>i 


Jln—k 


Then 

£> n+1 (x) = D(D n (x)) 

= x^ t(n, k)y k+l z 2n ~ k+1 +x^ kt(n, k)y k z 2n ~ k+2 + x J ^(2 n - fc)t(n, % fc+2 z 

fc>i fc>i fc>i 

Hence 

t(n + 1, k) = kt(n, k) + t(n , k — 1) + (2n — k + 2)f(n, k — 2). (9) 

By comparing Q with ([2]), we see that the numbers t(n, fc) satisfy the same recurrence relation 
and initial conditions as T(n,k ), so they agree. The assertion for D n (x 2 ) can be proved in a 
similar way. 

It follows from Leibniz’s formula that 


□ 


Hence 


D n (nv) = p 

D n (x 2 ) = jj C n ^jD k (x)D n ~ k (x), 

i —n \ / 


D n+ \x) = D n (xyz) = jh C‘f\D k (x)D n ~ k ( y z ) = £ D k (xy)D n ~ k (z). 

k =0 ' ’ k =0 ' ' 

Therefore, we can use Proposition [4] and Theorem [5] to get several convolution identities. 


ENUMERATION OF A DUAL SET OF STIRLING PERMUTATIONS BY THEIR ALTERNATING RUNS 7 


Corollary 6 . For n > 1, we have 

2(1 + x) n ~ l A n (x) = V (Er fc (x)T n _ fc (x), 
k =o ^ ' 

T n+ 1 (x) = x £ (f)T t (x)Bn-k(x 2 ), 

k = 0 ^ ' 

T n+ i(x) = Cf) S k (x)N n - k (x 2 ). 

k =0 ' ' 


( 10 ) 


Let T(x,z) = Yl^Lo R eca ll that the exponential generating function for A n (x) is 

given as follows (see [221 A008292]): 


. , . y - ^ , , t n 1 — X 

A(x,t) = ) A n (x)— = -77--r. 

v y ^ v 'n! l-xehi-^) 
n>0 


Combining (flOl) and (fill] . we get 


T(x,z) = 


e z(x-l){x+l) _j_ 


1 — x 2 


1 + X V e 2 z(z-l)0r+l) _ x 2 • 


( 11 ) 


( 12 ) 


From ([ 6 ]), we have 

E M «( 

Note that 


x 


n> o 


n! 


n> 0 

e z(a;-l)(a:+l) _|_ ,-j. 


1 \ Z n 


r 


1 — a; 2 


2/^| W g2 x(a:—l)(a:+l) _ ^,2' 


1 + X 

Therefore, from (|12lh we obtain 


i + E( x “ 1 ) n ( x+1 ) n_1 


n> 1 


71! 


n-1 , 

r n (s) = M n (x 2 ) + E UjM fc (x 2 )(x- l) n_fc (x + l) n “ fc_1 for n> 1. 

k =0 ' ' 


4. Concluding remarks 

In this paper, we show that the polynomial T n (x) has many similar properties to R n (x). In 
fact, there are more similar properties deserve to be studied. From the relation @ and the fact 
that A n (x) have only real zeros, Wilf [23j proved that R n {x) have only real zeros for n > 2. 
Moreover, it follows from a that all zeros of R n (x) belong to [—1,0], and the zeros of R n (x) 
separate that of R n+ \(x) (see |T7j Corollary 8 ]). 

Let /(x) and F(x) be two polynomials with only real coefficients. Suppose that /(x) and 
F(x) both have only imaginary zeros. We say that f(x) separates F(x) if degF = deg / + 2 
and the sequences of real and imaginary parts of the zeros of /(x) respectively separate that of 
F(x). In other words, let /(x) = a n"=i ( x + Pj + Qji)(x + Pj ~ Qji), an d let F{x) = hW"j =1 {x + 
Sj + tji)(x + Sj — tj i), where a,b are respectively leading coefficients of /(x) and F(x), p\ > p 2 > 

• • • > Pn-i, Qi > 92 > • • • > Qn- 1 ) Si > S 2 > • • • > Sn and t± > t 2 > ■ ■ ■ > t n . Then we have 

Sl>Pl> S 2 >P2> ■■■ > S n -1 > p n -l > S n , 

t\ > qi > t 2 > q 2 > ■ ■ ■ > tn -1 > Qn -1 > tn- 
Based on empirical evidence, we propose the following conjecture. 

Conjecture 7. Forn > 2, all zeros ofT n {x)/x are imaginary andT n (x)/x separates T n+ i(x)/x. 
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